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Abstract 

We study a class of stochastic evolution equations with a dissipative forcing 
nonlinearity and additive noise. The noise is assumed to satisfy rather general 
assumptions about the form of the covariance function; our framework covers 
examples of Gaussian processes, like fractional and bifractional Brownian 
motion and also non Gaussian examples like the Hermite process. We give 
an application of our results to the study of the stochastic version of a common 
model of potential spread in a dendritic tree. Our investigation is specially 
motivated by possibility to introduce long-range dependence in time of the 
stochastic perturbation. 
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1. Introduction 

In recent years, existence, uniqueness and further properties of solutions to stochas- 
tic equations in Hilbert spaces under dissipativity assumptions has been widely dis- 
cussed in the literature since similar equations play an important role in stochastic 
models of population biology, physics and mathematical finance (among others) , com- 
pare the monograph [5] for a thorough discussion. 

In this paper, by using semigroup methods, we shall discuss existence and uniqueness 
of mild solutions to a class of stochastic evolution equations driven by a stochastic 
process X which is not necessarily Gaussian. In the Hilbert space X we consider the 
following equation 

du(t) = Au(t) + F(u(t))dt + dX(t) ;|| . 
u(0) = Uq, 

where A and F satisfy some dissipativity condition on X and X is a general X- valued 
process that satisfies some specific condition on the covariance operator. 

Problems of the form of Equation (11. ip arise in the modeling of certain problems 
in neurobiology. In particular, in Section [5] we shall analyze a model of diffusion for 
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electric activity in a neuronal network, recently introduced in [3], driven by a stochastic 
term that is not white in time and space. Notice further that, motivated by this model, 
we are concerned with assumptions on the drift term which are not covered by those 
in 0. 

Assumption 1.1. The operator A : D(A) C X — » X is associated with a form 
(a, V) that is densely defined, coercive and continuous; by standard theory of Dirichlet 
forms, compare \19^ , the operator A generates a strongly continuous, analytic semigroup 
(S(t)) t >o on the Hilbert space X that is uniformly exponentially stable: there exist 
M > land lo > such that \\S(t)\\ L(x) < Me'" 1 for all t > 0. 

In the application of Section^ the operator A is not self-adjoint, as the correspond- 
ing form a is not symmetric; also, since V is not compactly embedded in X, it is easily 
seen that the semigroup generated by A is not compact hence it is not Hilbert- Schmidt. 

Assumption 1.2. F is an m-dissipative mapping with V C D(F) and F : V — ► X is 

continuous with polynomial growth. 

Let us introduce the class of noises that we are concerned with. We define the mean 
of a X valued process (-Xt)tefo T ] by mx ■ [0, T] — » X, mx(t) — EX t and the covariance 
C x : [0,T] 2 ^i!(X) by 

(C x (t, s)u, v) x = E \{X t - m x (t)v) x {X s - m x (s),u) x ] 

for every s,t£ [0, T] and for every u, v G X. 

Let Q be a nuclear self-adjoint operator on X (Q G £i(X) and Q — Q* > 0). It is 
well-known that Q admits a sequence (Xj)j>i of eigenvalues such that < Xj j and 
J2j ^j>i < oo- Moreover, the eigenvectors (ej)j>i of Q form an orthonormal basis of 
X. 

Let (x(t)) t £\o,T] be a centered square integrable one-dimensional process with a given 
covariance R. We define its infinite dimensional counterpart by 

Xt = Y,^ J {t)e ] t£[0,T], 

3=1 

where Xj are independent copies of x. It is trivial to see that the above series is 
convergent in L 2 (il; X) for every fixed t G [0, T] and 

E||X t |||= (TrQ)R(t,t). 

Remark 1.1. The process X is a X-valued centered process with covariance R(t, s)Q. 

Assumption 1.3. We will assume that the covariance of the process X satisfies the 
following condition: 

We will treat several examples of stochastic processes that satisfy Q1.2|) . The first 
two examples are Gaussian processes (fractional and bifractional Brownian motion) 
while the third example is non-Gaussian (the Hermite process). 
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Example 1.2. The process X is a fractional Brownian motion (fBm) with Hurst 
parameter H > |. We recall that its covariance equals, for every s,t £ [0, T] 

R(s,t) = ±( S 2H + t™-\s-t\ 2H ). 

In this case = 2H(2H — l)\t — s\ 2H 2 in the sense of distributions. Since R 

vanishes on the axes, we have for every s, t £ [0, T] 

R(s,t)= / dsi / ds 2 ^ — . 

Example 1.3. X is a bifractional Brownian motion with H £ (0,1), if £ (0,1] and 
2HK > 1. Recall the the bifractional Brownian motion (B t ' )t£[o,T] is a centered 
Gaussian process, starting from zero, with covariance 

R H '«(t, s) := R(t, s) = ± ((t™ + s 2H f \t s\ 2HK ) (1.3) 

with H £ (0, 1) and K £ (0, 1]. Note that, if K = 1 then B H1 is a fractional Brownian 
motion with Hurst parameter H £ (0, 1). 

Example 1.4. A non Gaussian example: the Hermite process: The driving process is 
now a Hermite process with selsimilarity order H £ (^,1). This process appears as a 
limit in the so-called Non Central Limit Theorem (see [7] or [25] ). 

We will denote by {Z\ [ q ' )te[o,i] the Hermite process with self- similarity parameter 
H £ (1/2, 1). Here q > 1 is an integer. The Hermite process can be defined in two 
ways: as a multiple integral with respect to the standard Wiener process (Wt)te[o,i]y or 
as a multiple integral with respect to a fractional Brownian motion with suitable Hurst 
parameter. We adopt the first approach throughout the paper: compare Definition ^. 21 
below. 

In Section [3] we treat the stochastic convolution process 

W A (t) = f S(t-s)dX s . (1.4) 
Jo 

It is the weak solution of the linear stochastic evolution equation dY(t) = AY(t) dt + 
dX(t). Our aim is to prove that is a well-defined X- valued, mean square continuous, 
^-"t-adapted process. We strenghten Assumption (|1.2[) by imposing the following. 

Assumption 1.4. Let X be given in the form 

where Xj, ej and Xj(t) have been defined above. Suppose that the covariance R of the 
process (Xt)te[o,T] satisfies the following condition: 

< Cl \t-s\ 2H - 2 +g(s,t) 

for every s,t £ [0,T] where \g(s,t)\ < C2(st) 13 with [3 £ (—1,0), H £ (4,1) and ci,C2 
are strictly positive constant. 



OR 



dsdt 
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Remark 1.5. The fractional Brownian motion and the Hermite process satisfy As- 
sumption [TJ2] with g identically zero. In the case of the bifractional Brownian motion 
the second derivative of the covariance can be divided into two parts. Indeed 

g(u, v) =a\u- v\ 2HK - 2 + c 2 {u 2H + v 2H ) K - 2 (H 2ff ~ 1 := «) + v). 

The part containing gx can be treated similarly to the case of the fractional Brownian 
motion. For the second term, note that 

u 2H + V 2H > 2{uv) H and {u 2H + y 2H )K -2 < 2 K ~ 2 (uv)"^. 

So, 

| 52 (w,w)| < cst.{uv) HK - 1 . 

In conclusion, Assumption II .41 is satisfied with (3 — HK — 1 S (—1, 0). 

Our first main result is the following theorem concerning the regularity of the 
stochastic convolution process under the Assumptions 11.11 and 11.41 

Theorem 1.5. In the above framework, fix a G (0,11). Let W& be given by \l-4-\ ). 
Then Wa exists in L 2 ([0,T] x f2;X) and it is Tt-adapted. 
For every 7 < a and e < a — 7 it holds that 

W A eC a ^- £ ([0,T};D(-Ay); 

in particular for any fixed t E [0, T] the random variable W&(t) belongs to D(-A)' 1 . 

Now we consider the solution of the stochastic evolution equation (|1.1[) . Wc consider 
generalized mild solutions in the sense of [SJ Section 5.5]: an X- valued continuous and 
adapted process u — {ut, t > 0} is a mild solution of (II. ip if it satisfies Pr-a.s. the 
integral equation 

u(t) = S(t)uo+ f S{t-s)F{u{s))ds + W A {t). (1.5) 
Jo 

Theorem 1.6. In our setting, let uq € D(F) (resp. Uq G X). Then there exists a 
unique mild (resp. generalized) solution 

u e L^Cfl; C([0, T]; X)) n h 2 ^- L 2 ([0, T]; V)) 

to equation which depends continuously on the initial condition: 

E\u(t;u ) - u(t\Ui)\ x < C\u - ui| x . (1.6) 

Remark 1.6. Even for a Wiener perturbation, this result is not contained in the 
existing literature since we does not assume any dissipativity or generation property 
of A on V, compare \S\ Hypothesis 5.4 and 5.6]. 

With this result at hand, we can solve the model of a complete neuronal network 
recently proposed in It is well known that any single neuron can be schematized 
as a collection of a dentritic tree that ends into a soma at one end of an axon, hence 
as a tree in the precise sense defined within the mathematical field of graph theory. By 
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introducing stochastic terms we can model the chaotic fluctuations of synaptic activity 
and post-synaptic elaboration of electronic potential. There is sufficient sperimental 
evidence that, in order to capture the actual behaviour of the neurobiological tissues, 
infinite dimensional, stochastic, nonlinear reaction-diffusion models are needed. 

Previous models used simplified version of the neuronal network or just concentrate 
on single parts of the cell: compare [T] for a thorough analysis of the FitzHugh Nagumo 
system on a neuronal axon or [5] for the analysis of the (passive) electric propagation 
in a dendritic tree in the subtreshold regime. In our model, instead, we are based based 
on the deterministic description of the whole neuronal network that has been recently 
introduced in Q; therefore, we avoid to sacrify the biological realism of the neuronal 
model and, further, we add a manifold of different possible stochastic perturbations 
that can be chosen as a model for the enviromental influence on the system. Notice 
that already in [5] a fractional Brownian motion with Hurst parameter H > i was 
chosen in order to model the (apparently chaotical) perturbance acting on a neuronal 
network. This choice is not a premiere in neuroscience, since different considerations 
show that real inputs may exhibit long-range dependence and self-similarity: see for 
instance the contributions in [22j Part II]. 

2. Wiener Integrals with respect to Hilbert valued Gaussian and 
non-Gaussian processes with covariance structure measure 

In this section we discuss the construction of a stochastic integral with respect to 
the process (Xt)t£[o,T], which is not necessarily Gaussian. 

Since our non-Gaussian examples will be given by stochastic processes that can be 
expressed as multiple Wiener-Ito integrals, we need to briefly recall the basic facts 
related to their constructed and their basic properties. 

2.1. Multiple stochastic integrals 

Let (Wt) te [ 0j T] be a classical Wiener process on a standard Wiener space (fi,.? 7 , P). 
If / 6 L 2 ([0, T] n ) with n > 1 integer, we introduce the multiple Wiener-Ito integral of 
/ with respect to W . The basic reference is the monograph [18]. Let / £ S m be an 
elementary function with m variables that can be written as 



where the coefficients satisfy c%i,...i m = if two indices ik and %\ are equal and the sets 
Ai e B([0, T]) are disjoints. For such a step function / we define 



where we put VF([a, b]) — Wb — W a - It can be seen that the application I m constructed 
above from S m to L 2 (fl) is an isometry on «S TO , i.e. 




l\ ,. . . ,2 



im(/)= c lu ..., im W(A n )...W(A z J 



^[In(f)I m {g)] = n\{f,g) L 2 {[0T]n) if m = n 



(2-1) 



and 



E[I n (/)J m ( ff )] = Oifm^ra. 
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Since the set S n is dense in L 2 ([0,T] n ) for every n > 1, the mapping /„ can be 
extended to an isometry from L 2 ([0, T] n ) to L 2 (il) and the above properties (|2.1D hold 
true for this extension. 

We recall the following hypercontractivity property for the L p norm of a multiple 
stochastic integral (see [131 Theorem 4.1]) 

E|/ m (/)| 2m <c m (E/ m (/) 2 )™ (2.2) 

where c m is an explicit positive constant and / G L 2 ([Q,T] m ). 

2.2. Wiener integrals: the one-dimensional case 

The idea to define Wiener integrals with respect to a centered Gaussian (or non 
Gaussian) process (Xt)te[o,T] is natural and standard. Denote by R(t,s) = E(X t X s ) 
the covariance of the process X. Consider £ the set of step functions on [0, T] defined 

as 

n-1 
i=0 

where n : = to < t\ < . . . < t n = T denotes a partition of [0, T] and Cj are real 
numbers. For a such / it is standard to define 

n-1 

/ (/) = E c ^ {x u+1 -x u ). 

It holds that 

n-1 

i,j=0 
n-1 

= E CiCj (R(ti + i,tj + i) — R(t i+ i,tj) — R(ti,tj + i) + R(U,tj)) . 

i,j=0 

The next step is to extend, by density, the application I : £ — * L 2 (il) to a bigger space, 
using the fact that it is an isometry. This construction has been done in [TU] and we 
will describe here the main ideas. In particular, we shall see that the construction 
depends on the covariance structure of the process X\ the covariance of X should 
define a measure on the Borel sets of [0, T) 2 . The function R defines naturally a finite 
additive measure \i on the algebra of finite disjoint rectangles included in [0,T] 2 by 

fi(A) = R(b, d) + R(a, c) - R(a, d) - i?(c, b) 

if A = [a,b) x [c,d). 

In order to extend the Wiener integral to more general processes, we assume that 
the covariance of the process X satisfies the following condition: 

(M) ^S eLl([o ' T]2) - (2 - 4) 

(compare with Assumption |1.3|l . This is a particular case of the situation considered 
in [lOj where the integrator is assumed to have a covariance structure measure in the 
sense that the covariance R defines a measure on [0, T] 2 . 
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We have already seen some examples of stochastic processes that satisfy (|1.2p : the 
fractional Brownian motion with Hurst index bigger than i , the bifractional Brownian 
motion with 2HK > 1 and the Hermite process, for instance. 

The next step is to extend the definition of the Wiener integral to a bigger class of 
integrands. We introduce \Tt\ the set of measurable functions / : [0, T] ~ > R such that 



T P T 



\f(u)f(v) 



d 2 R 



dudv 

On the set \TL\ we define the inner product 



(u,v) 



du dv < oo. (2-5) 



(J.h)u= / [ f(u)h(v)-^(u,v)dudv (2.6) 



JO 



and its associated seminorm 



d 2 R 

f( u )f( v ) g u Q ( u -> v ) du dv - 



We also define 



\\H\ 



T 



d 2 R 



dudv 



(u,v) 



dudv. (2.7) 



It holds that £ C \H\ and for every f,h € £ 

EI(ff =E\\f\\n- (2-8) 
The following result can be found in [TP] , 

Proposition 2.1. The set £ is dense in \H.\ with respect to |H|i^i and in particular to 
the seminorm The linear application $ : £ — ► L 2 (fl) defined by 

tp — > I(ip) 

can be continuously extended to \7i\ equipped with the ||-|| w -norm. Moreover we still 
have identity \2.8\) for any tp 6 \H\. 

We will set tpdX = $(tp) and it will be called the Wiener integral of tp with 
respect to X. 

We remark below that if the integrator process is a process in the nth Wiener chaos 
then the Wiener integral with respect to X is again an element of the nth Wiener 
chaos. 

Remark 2.1. Suppose that the process X can be written as X t = Ik{L t {-)) with k > 1 
and L t e L 2 ([0,T] k ) for every t E [0,T]. Then for every tp 6 |W| the Wiener integral 
J Q T y> dX is also in the fcth Wiener chaos. Indeed, for simple functions of the form (|2.3|) 
it is obvious and then we use the fact that the kth Wiener chaos is stable with respect 
to the L 2 convergence, that is, a sequence of random variables in the fcth Wiener chaos 
convergent in L 2 has as limit a random variable in the fcth Wiener chaos. 
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Remark 2.2. Assumption 11.41 implies condition (|1.2|) . In particular the process X 
whose covariance satisfies Assumption 11.41 have a covariance structure measure and 
the Wiener integral ipdX exists for every ip 6 \H\. Indeed, 

! ! ^rS-(M) dsdt< Cl [ [ \s-t\ 2H - 2 dsdt + c 2 [ [ (stfdsdt 
Jo Jo osdt J a J a J a J 

< c ( T 2H + T 2(/3+i; 



Let us discuss now some examples. Firstly we refer to Gaussian processes (fractional 
and bifractional Brownian motion). 

Example 2.3. The case of the fractional Brownian motion with H > ^- In this case 
\H\ is the space of measurable functions / : [0, T] — > R such that 

T 

\f(u)f(v)\u - v\ 2H ~ 2 dudv < co. 

On the other hand, for this integrator one can consider bigger classer of Wiener 
integrands. The natural space for the definition of the Wiener integral with respect to 
a Hermite process is the space TL which is the closure of £ with respect to the scalar 
product 

(l[o,t]j l[o,s])« = s )- 

We recall that TL can be expressed using fractional integrals and it may contain 
distributions. Recall also that 

L 2 ([0,T}) c ([0,2*1) c \H\ c H 

The Wiener integral with respect to Hermite processes can be also written as a Wiener 
integral with respect to the standard Brownian motion through a transfer operator 
(see e.g. [B]). 



Example 2.4. The bifractional Brownian motion with 2HK > 1. Recall the the 
bifractional Brownian motion [B^ ,K ) t e[o.T] is a centered Gaussian process, starting 
from zero, with covariance 

R H ' K (t, s) := R(t, s) = -L ({t 2H + s 2H f \t s\ 2HK ) (2.9) 

with H e (0, 1) and K G (0, 1]. 

We can write the covariance function as 

R{s\,s 2 ) = Ri(si,s 2 ) + i?2(si, S2), 

where 

Ri(i,n) - ^ {4 H + 4 H ) K ~ (4 HK + si HK ) 

and R 2 (s u s 2 ) = -i| S2 - Sl \ 2HK + s\ HK + s% HK . 
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Wc therefore have 

d 2 Ri _ 4H*K(K - 1) , 2H 2H ,k-2 2H _ X 2H _, 

d^d7 2 ~ 2^ ^ + * 2 > Sl ° 2 ■ 

Since R\ is of class C 2 ((0,T] 2 ) and q s qI is always negative, R\ is the distribution 

function of a negative absolutely continuous finite measure, having g s for density. 

Concerning the term i? 2 we suppose 2HK > 1. The part denoted by R2 is (up to 
a constant) also the covariance function of a fractional Brownian motion of index HK 
and = 2HK(2HK - 1) \s x - s 2 \ 2HK ~ 2 which belongs of course to Z^QO.T] 2 ). 

We also recall that the bifractional Brownian motion is a self-similar process with self- 
similarity index HK, it has not stationary increments, it is not Markovian and not a 
semimartingale for 2HK > 1. 

A significant subspace included in \H\ is the set £ 2 ([0,T]); if K = 1 and H = 5, 
there is even equality, since X is a classical Brownian motion (see [10] V 

Let us now give a non-Gaussian example. 

Example 2.5. TTie Hermite process Z ( * q,H " > := Z of order q with self similarity order 
H. 

The fractional Brownian process (Bf )tg[o,i] with Hurst parameter ii G (0, 1) can 
be written as 

Bf - / X H (M)dW/, t g [0,1] 
Jo 

where (Wt, t G [0, T]) is a standard Wiener process, the kernel K H [t, s) has the expres- 
sion cms 1 ' 2 -" $l{u - s)"-^ 2 ^- 1 ' 2 du where t > s and c H = ( m "?H,H%2) ) ^ 
and /?(-, •) is the Beta function. For t > s, the kernel's derivative is 9 ^ (t, s) = 

ch (f) 1 ^ 2 H (t ~ s) H ~ 3 / 2 . Fortunately we will not need to use these expressions 
explicitly, since they will be involved below only in integrals whose expressions are 
known. 

We will denote by (Z ( t q - H) ) te[QA] the Hermite process with self-similarity parameter 
H G (1/2,1). Here q > 1 is an integer. Let us state the formal definition of this 
process. 

Definition 2.2. The Hermite process (Z^ q ' )te\o,i] °f order q > 1 and with self- 
similarity parameter H G (5, 1) is given by 

Z^ H) = d{H) f ... f dW yi . . . dW yq ( /* d x K H ' (u,y x )... d 1 K H ' («, „,)d« ) 

JO JO \J yi V...Vy q J 

te[o,i] (2.10) 

where isT^ is the usual kernel of the fractional Brownian motion and 

TT 1 

H' = l + <^ (2H' - 2)q = 2H - 2. (2.11) 
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Of fundamental importance is the fact that the covariance of Z^ H ^ is identical to 
that of fBm, namely 



E 



Z 



(q,H) 7 (q-.H) 



1 



2H | 2H 



\t-»ry 



The constant d(H) is chosen to have the variance equal to 1. We stress that Z^ q,H ^ is 
far from Gaussian for q > 1, since it is formed of multiple Wiener integrals of order q 
(see also [26 ). 

The basic properties of the Hermite process are listed below: 

• the Hermite process Z^ is ii-self-similar and it has stationary increments. 

• the mean square of the increment is given by 



z(q,H) _ 7 ( q ,H) 



= t-s 



2H. 



(2.12) 



as a consequence, it follows will little extra effort from Kolmogorov's continuity 
criterion that Z^ q ' H ^ has Holder-continuous paths of any exponent S < H. 
• it exhibits long-range dependence in the sense that 



J2^[z{ q ' H) (z^-z^) 



n>l 

In fact, the summand in this series is of order n 2H ~ 2 . This property is identical 
to that of fBm since the processes share the same covariance structure, and the 
property is well-known for fBm with H > 1/2. 

• for q = 1, Z^-' H ^ is standard fBm with Hurst parameter H, while for q > 2 
the Hermite process is not Gaussian. In the case q = 2 this stochastic process is 
known as the Rosenblatt process. 

In this case the class of integrands Ti. is the same as in the case of the fractional 
Brownian motion. We will also note that, from Remark 12.11 the Wiener integral 
with respect to the Hermite process tpdZ is an element of the fcth Wiener chaos. 
Moreover, it has been proven in [12] that for every ip <E \H\ we have 



f(u)dZ(u) 



f ... [ I(f)(y 1 ,y 2 ,... : y k )dB(y 1 )dB(y 2 )...dB(y k 
Jo Jo 



where (-B*)*6[o,t] is a Wiener process and we denoted by / the following transfer 
operator 



I(f)(y u y 2 ,...,y k )= I f(u)d 1 K H '(u,y 1 )...d 1 K"(u,y k )du 
JyiV-VVk 



where H' is defined by (|2.1ip . 
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2.3. The infinite-dimensional case 

Let 

•v.. ^r X A> /: /: (/ *G[o,n 

be a X-valued centered process with covariance i?(t, s)Q. 

Let G : [0, T] — > L(X) and let (e 3 -)j>i be a complete orthonormal system in X. 
Assume that for every j > 1 the function G(-)e 3 belongs to the space \H\. The we 
define the Wiener integral of G with respect to X by 



[ T GdX = Y,^ [ T G(s)e j dx j (s) 
Jo 7 >i ^0 



where the Wiener integral with respect to dxj has been defined above in paragraph 



Remark 2.6. The above integral is well-defined as an element of L 2 (f2;X) and we 
have the bound 









E 


/ GdB 


< Tr(Q) j T j 




Jo 


Jo Jo 



d 2 R 



dudv 



(u,v) 



du dv 



< Tr(Q) \\\\G(-)\\ L (v)\\\ n \ • 

3. The stochastic convolution process 

There exists a well established theory on stochastic evolution equations in infinite 
dimensional spaces, see Da Prato and Zabcyck [S], that we shall apply in order to show 
that Eq. (|1.5[) admits a unique solution. Let us recall from Assumption 11.11 that A is 
the infinitesimal generator of a strongly continuous semigroup (S(t)) t >a, on X that is 
exponentially stable. 

In this setting, we are concerned with the so-called stochastic convolution process 



w A (t) = f s(t-s)dx s . 03) 

Jo 

It is the weak solution of the linear stochastic evolution equation dY(t) = AY(t) dt + 
dX(t). Our aim is to prove that is a well-defined mean square continuous, J^-adapted 
process. Let us make the following assumption: 

Proposition 3.1. Assume that the covariance function R satisfies ll.ty) . Then, for 
every t G [0, T], the stochastic convolution given by \l-4\ l exists in L 2 ([0,T];X) and it 
is Tt adapted. 

Proof. We have that, by using the exponential stability of the semigroup S(t) (see 
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Assumption II. 1|) 



E 



S(t-s)dX(s) 



<Tr(Q) / / \\ S (t-n)\\ HX) \\S(t-v)\\ L{x) 
Jo Jo 

d 2 R 



d 2 R 



dudv 



(u, v) 



du dv 



T r T 



< M 2 Tr(Q) f f 
Jo Jo 



g -ui(t-ii) e -u(t-v) 



dudv 



(u,v) 



dudv < oo. 



The fact that W& is adapted is obvious. 



The next step is to study the regularity (temporal and spatial) of the stochastic 
convolution process. This will lead to a study of an infinite sum of random variables 
with independent but not necessarily Gaussian summands (they are elements in a fixed 
order Wiener chaos). Let us recall the following result from [TTJ Theorem 3.5.1, page 
76 and Theorem 2.2.1, page 32]. 

Proposition 3.2. Let X be a Hilbert space. 

a) Let p > 4 and Xi, . . . , X n be zero mean, independent X valued random variables. 
Then 



E 



< c„ 



2\ 2 



+ (E||X„|||v(E||X n _ 1 |||v(...V]E||X 1 |||)))5 



b) Let p > and X\, X2, ■ ■ • , X n , ... be a sequence of independent X valued random 

variables. If the series ^ Xi converges almost surely to a random variable S and 

i>i 

for some t > 



^E||X i |||l ( || X< || x>t) <oo 



(3.1) 



then E \\S\\Z < 00 and E \\S n — S\\Z —> where S n — ^2 Xi 



Remark 3.1. It is not difficult to see that the point a) above implies that 



E 



E E ii^n p 



(3.2) 



The following lemma is the main tool to get the regularity of the stochastic convo- 
lution process Wa- 



Dissipative stochastic evolution equations driven by general noise 



13 



Lemma 3.3. Let (^t)te[o,T) a stochastic process whose covariance R satisfies Assump- 
tion \1.4\ Denote, for every a E (0, 1), 

Y a (t)= [ {t-u)- a S(t-u)dX u , te[0,T}. 



Then for every a E (0, H), Y a belongs to LP ([0, T]; X). 

Proof. Suppose first that X is Gaussian. Then, in order to show that Y a is in 
LP ([0, T];X) it suffices to proves that it is in L 2 ([0, T]; X). We have 



\\\Y a (t)\\i=K 



^ C E A * f f{t-u)- a {t-vr a \\s{t 

• >x Jo Jo 



u)ej\\ x \\S(t - v)ej\\ 
d 2 R 



<C(TtQ) / / {t-uy^it-vy^e-^-^e-^-^ 



o Jo 

t r t 



dudv 
d 2 R 



dudv 



(u,v) 
(u,v) 



du dv 



du dv 



<C(TrQ) / / (t - u)- a (t - v^e-^-^e-^^^u - v\ 2H ~ 2 dudv 



Jo 

+ C(TrQ) f [ (t-u)- a {t~v)- a e- Ul{t - u) e-^ t - v \uvf dudv 
JO Jo 

■=h + h. 

Concerning the term 2i, we can write 

h<2C(TrQ) f [ {t-u)- a {t-vy a \u-v\ 2R - 2 dvdu 
Jo Jo 



<C(TrQ) fu-^u 211 - 1 [ z~ a {\~ z) 2H - 2 dzdu^C{YiQ) f 
Jo Jo Jo 



u- 2a u 2H ~ l 



du 



where we used the change of variable — = z. The last quantity is clearly finite if and 
only if a < H . Concerning Li we have 



h < (J {t-uye-^-^u du^j 



and this is always bounded by a constant (depending only on T) using the hypothesis 
imposed on a and (3. We obtain thus the bound 

E\\Y a (t)\\ 2 x <C = C T 



for every t E [0,T]. 
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Let us assume now that X is not Gaussian and it belongs to the fc-th Wiener chaos 
with k > 2. The process X can be written as 

x t = v'Vi^w 

where Xj is an element of the k th Wiener chaos with respect to the Wiener process 
Wj and (wj)j>i are independent real one-dimensional Wiener processes. Then 



i>i Jo 



is also an element in the k th Wiener chaos (in the sense that every summand is in the 
k th Wiener chaos with respect to Wj). Note also that, using the above computations 
from the Gaussian case we obtain 

E\\Y a (t)\\l<C = C T 

for every t £ [0, T]. Denote by 

n „t n 

Sn,Y a (t) - E V^' / (t - u)- a S(t - u) dx 3 := £ ^/aJAj. 
By Proposition 13. 2[ point a) and Remark l3~Tl we have 







n 


2\ 5 


(Eii^jn* < 


E 










3=1 


X 



E 



E 



Using the hypercontractivity property of multiple stochastic integrals ()2.2|1 . we get for 
every i = 1, .., n 



A,; A 



= Aj 



(< - u)- a S(t - u) d Xj 



(t - u)- a S{t - u) dxj 



2\ 2 



< Cp, T A? . 



As a consequence, since < A, j 

E||S„, K ,(f)|||<c p , T I |J> 



< c 



p,T< 



(3.3) 



Now, since for every t the sequence S n ,y Q is convergent in L 2 (f2;X) as n -* oo we can 
find a sequence which converges almost surely. This subsequence will be again denoted 
by S n ,Y a - By Proposition 13.21 point b), since 



UaTA l (VX r Ai>t) <^A?E||AJ|<c p , T 
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we obtain that for every t the random variable Y a (t) belongs to L p ([0, T]\ X) and 
E \\S rh Y a (t) — ^a(^)|lx ^n^oo 0. Letting now n — » oo in (|3.3p we obtain that 

Ep r e ,(i)||£<c p ,T. 

and this finishes the proof. 

Proposition 3.4. Suppose that X satisfies Assumption \ l.J\ and fix a £ (0,H). Let 
Wa be given by jl.4\ l- Then for every 7 < a and e < a — 7 if ZioZcfe i/iat 

^ A £C tt ^- £ ([0,r];I?((-AP)). 

In particular for any fixed t £ [0,T] i/ie random variable W^it) belongs to D((— A) 7 ). 
Proof. For a, 7 £ (0, f ), p > 1 and V £ £ p ([0, T], X) we define 

P Q ^(i) = sm ( a7r ) / (i - u) a -\-A)~<S(t - u)ip(u) du. 
Jo 

Then, if a > 7 + i it holds that 

R an £ L (Lf([0,r];X);C a - 7 -^([0,T];^((-A) 7 ))) 
It is standard to see that 

{-A)~>X(t) = {R an Y a ){t) 

where Y a (t) = f*(t-u)- a S(t-u) dX(u). Since by the above lemma Y a £ Z,*>([0,T];X) 
the conclusion follows. 

Next we will regard further properties of the stochastic convolution process (|1.4p . 
We are concerned with the LP norm of its supremum and with its regularity with 
respect to the time variable. In the Gaussian case the proofs basically follow the 
standard ideas from [6j Chapter 5], while in the non-Gaussian case, the results are new 
and they involve an analysis of the LP moments of multiple Wiener-Ito integrals. 

Lemma 3.5. Assume that {X t )te[Q,T] satisfies Assumption ] 1.4\ and let Wa be given by 
\l-4\l - For any p > jj, we have 

E sup \\W A (t)\\l<C. 
te[o,T] 

Proof. Note that for t £ [0, T] and < a < H 

W A (t) = ^^ ( Sit-s^t-sY^Z^ds 
t Jo 

with Z a (s) — J Q S S(s — u)(s — u)~ a dX u . By Holder's inequality with P > > jj 

E sup ||W A (t)|||<E / \\Z a (s)\\l ds. 
te[o,T] Jo 
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Now, since < a < H, 

2 <r f T -uj(t-u) -u(t-v) d2R 



E\\Z a (,)\&<C e-^-^e 



dudv 



(u,v) 



dudv < C 



and by using the computations in the proof of Lemma 13.31 and the hypercontractivity 

\bX 



property for multiple stochastic integrals (|2.2[) we get E ||^a:(s)|lbx — C ■ 



Let us now state our result concerning the regularity of W A with respect to the time 
variable. 

Proposition 3.6. Fix a G (0, H A ((3 + 1)). Then the process W&(-) has a Holder 
continuous paths. 

Proof. We will use Kolmogorov's continuity criterium for Hilbert valued stochastic 
processes (see Theorem 3.3]). To this end, we need the evaluate the increment 
W A {t) - W A (s). We can write 



W A (t) - W A (s) = J2 I Sit - u) ej d Xj (u) 

+ y2v^i (S(t-u) - S(s-u))ejdxj(u) 
„-^i Jo 



3>1 

■=h + h. 



Concerning the first term, we get from Assumption II. 41 

Kll<(TrQ)cxj j e-^-^e-^-^lu-vfx-idudv 
+ {TrQ)c 2 f f e'^-^e-^-^iuvf dudv 



~ C Us is \ U ~~ V \ 2H 2dudV + f s j g ( uv f dudv 

< C (\t - S \ 2H + \t - s\W+V). 

Following the proof of Theorem 5.1.3], we obtain Eif < C\t-s\ 2 "> for any 7 G (0, 1). 
As a consequence 

E ||W A (t) - Wa{s)\\1 < C (|t- s\ 2H + \t- s\W+V + \ t _ a |*r) 

and by (|2.2[) we will have that (as in the proof of Lemma |3~3"]) for every s close to t 

E ||W A (t) - W A (s)\\l <C p (\t- s\*> H + \t- s |^ +1 )) 

and this bound will imply the existence of an a-H61der continuous version of T^a- 

Remark 3.2. In the case of the fractional Brownian motion and of the Rosenblatt 
process the order of continuity is H. For the bifractional Brownian motion, since 
fj + 1 = HK, the stochastic convolution is HK Holder continuous. 
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4. Existence and uniqueness of the solution 

Let us first introduce the spaces where the solution will live. 

Definition 4.1. Let Lj?(ri; C([0, T]; X)) denote the Banach space of all ^-measurable, 
pathwisc continuous processes, taking values in X, endowed with the norm 



l x llz^(n ; C([o,T];X)) = I E t ^ T H X (*) 



\ 1/2 
|2 \ 



while L%(fl;L 2 ([0,T];V)) denotes the Banach space of all mappings X : [0,T] -> V 



such that is .Ft-measurable, endowed with the norm 



X \\l*(Q;L 2 (_[0,T]:V)-) = I E / ll^Wllyd* 



'i' 



1/2 



,2 







We are concerned with Eq. (jl.ip that we mean to solve in mi/d /orm: a process 
w e L^(fi;C([0,T];X)) n L%{tt; L 2 ([0, T]; V)) is a solution to Eq. (TTTJl if it satisfies 
P-a.s. the integral equation 

u(t) = S{t)u + [ S(t-a)F(u(a))da + W A (t), te[0,T}. 
Jo 

The strategy of the proof is classical, compare [SJ Theorem 5.5.8]: we consider the 
difference {u(t) — WA(t))te[o,T] and we prove that it satisfies the mild equation and it 
belongs to the relevant spaces. 

4.1. Existence of the solution for deterministic equations 

Let us consider the following evolution equation 

^y(t)=Ay(t) + F(z(t) + y(t)) 
2/(0) = u Q , 

where A and F satisfy the dissipativity condition on X stated in Assumptions 11.11 
and 11.21 and z is a trajectory of the stochastic convolution process, which satisfies the 
regularity conditions stated in Theorem 11.51 

z <z c a -T e {%T];D(-Af) . 

The construction in this section is based on the techniques of Section 5.5]; notice 
however that we are concerned with a different kind of stochastic convolution and we 
do not impose any dissipativity on the operators A and F on the space V. 

Remark 4.1. The key point in the following construction is the observation that 
V = D{{— A) 1 / 2 ), compare Remark 15.31 Further, in this case we impose the following 
bound: h < H . Therefore, we can and do assume that 



for arbitrary e > 0. 

Now, notice that the assumption on F implies that F : V — > X is continuous, hence 
the process (-F(z(0))te[o,T] is continuous and satisfies sup tg [ T ] ||F(z(f))||x < +oo. 
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Let us introduce the Yosida approximations F a of F. It is known that F a are 
Lipschitz continuous, dissipative mappings such that, for all it G V, it holds F a (u) — > 
F(u) in X, as a — > 0. 

In this part, we are concerned with the following approximation of Eq. (14. lft : 

^-y a (t) = Ay a (t) + F a {z[t) + y a (t)) 
y a (0) = u . 

Lemma 4.2. Let x G X. Then, for any a > there exists a unique mild solution 
y a (t,x) to Eq. J^.£| ) such that 

y a eC([0,T];X)nl 2 ([0,r];¥). 

Proof. Since F a are Lipschitz continuous, the existence of the solution to (14. 2|) is 
standard. It remains to prove the existence of a estimate that is uniform in a. 

By the assumptions on A there exists lu > such that (Am, u) < —u\\u\\y, compare 
also Remark 15. 4) using the dissipativity of F we have 

TrlbaWll! = 7rlK||l+ / (Ay a (s),y a (s)) x ds + [ (F a (z{s) + y a (s)),y a (s)) x ds 



< -^\\uo\\l - ^ I ||2/aO)Hv ds + / (F a (z(s)),y a (s))xd& 



<dl«o||x-w/ hMU^+T sup \\F(z(t))\\i+ / ||y a (fl)||£ds 
z Jo te[o,T] Jo 

which implies, by an application of Gronwall's lemma, that 

sup fj||»a(t)||| + w / ||i/ a («)||^da > ) <C(T,u 0> z). (4-3) 



te[o,T] \2 

Notice that the constant on the right-hand side is independent of a. 
Lemma 4.3. For every a > 0, Uq,Ui G X, Zio/cfe 

sup lll/^W-^WIIl^CHuo-Uilll. (4.4) 
te[o,T] 

Proof, us consider the difference y^ a {t) — y^if), for x,x £ H: 

^ [<C («) - z/. 1 (*)] = a [«) - (*)] + [*"«(*(*) + CW) - F„(*(t) + !#(«))] 

hence 

1 IjC (*) - ?C (*) I II = IK - III + 2 / (yl(yS° (*) ~ vT (*))> y u a ° (a) - (*)> ds 

Jo 

+ 2 / (^(^(s))-^^ 1 ^)),^ ^)-^ 1 ^)) ds 
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and therefore 

111/? (<) ~ C Will < IK - III - 2c / («) - j# (s)\\l da. 

Jo 

Applying Gronwall's lemma we obtain 

IICW-^WIIx^e-^IK-^illx- (4-5) 
Lemma 4.4. The sequence (y a )a>o is a Cauchy sequence in C([0, T]; X)ni 2 ([0, T]; V). 
Proof. Let a, /? > 0. Then we compute 

^ M*) - y/j(t)] = A [y a {t) - yp(t)] + [F a (z(t) + y a (t)) - F p (z(t) + y (t))} 
Now, let us recall that 

(F a (x) - F p { V ) lX - y) x <(a + 0) [\F a (x)\ 2 + \F p {y)\ 2 } 
for all i,i/eV,a,(?>0 (compare Proposition 5.5.4]; it follows that 

1 f* 

<(« + /?) /Va(2(s)+y Q (.s))||| + ||^(z( S )+yMs))llx ds - ( 4 - 6 ) 
Jo 

Since F : V — > X is continuous, it follows that for some contant L > 



||F a (^)+2/a(s))||x< ||f(*(«)+lfa(«))llx 

<L ||z(.s) +2/Q ( s )||2 <2L [||z(a)|ft+lly«0OII$] 

hence by using estimate (I4.3|) 

||F Q (z( s ) + y a (s))|£ + \\F p (z{s) + yp(s))\\ 2 x ds 

< 2LT||z||p ([o r] . v) + C(T, u , 2, w, L) 



o 



is bounded by a constant that does not depend on a and /3. If we put the above 
estimate in (14.61) we obtain 



- sup \\y a (t)-yp(t)\\ x + w \\y a (s)-yp(s)\\ x ds<C(a + P) 
z te[o,T] Jo 

which easily implies the thesis. 

Theorem 4.5. For any z G C([0,T];V) i/iere exists a unique solution (y(t))te[0T] to 
Eq. 

2/GC([0,T];V)nL 2 ([0,T];X) 
and it depends continuously on the initial condition Uq G X. 
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Proof. Since y a is a Cauchy sequence in C([0, T];V) D L 2 ([0,T];X) it converges to 
a unique function y in the same space; it remains to show that (y(i))tg[o,T] actually 
solves (|4.1[) . Also, the continuous dependence on the initial condition follows from the 
same property proved for the approximating functions y a , since the estimate in (|4.4I) 
does not depend on a and it is conserved at the limit. 

By the claimed convergence of y a , since J a is a sequence of continuous mapping that 
converges to the identity, it holds that J a (y a (s)) — » y(s) € V a.s. on [0,T]. Therefore, 
by the continuity of F, it follows that 

F a {z(s) + y a {s)) ^ F(z(s) + y(s)) eX a.s. on [Q,T]. 

Now we use Vitali's theorem (the Uniform Integrability Convergence Theorem, com- 
pare Theorem 9.1.6]), to conclude that 

S{t-s)F a {z(s)+y a {s))ds^ / S{t-s)F(z{s)+y{s))ds. 

Jo 

5. A network model for a neuronal cell 

In this paper we aim to investigate a mathematical model of a complete neuron which 
is subject to stochastic perturbations; for a complete introduction to the biological 
motivations, see [9]. Our model is based on the deterministic one for the whole neuronal 
network that has been recently introduced in [4]; we shall borrow from this paper the 
basic analytical framework for the well-posedness of the problem. 

We treat the neuron as a simple graph with different kind of (stochastic) evolutions 
on the edges and dynamic Kirchhoff-type condition on the central node (the soma). 
This approach is made possible by the recent development of techniques of network 
evolution equations; hence, as opposite to most of the papers in the literature, which 
concentrate on some parts of the neuron, could it be the dendritic network, the soma 
or the axon, we take into account the complete cell. 

In this paper, we schematize a neuron as a network by considering 

• a FitzHugh-Nagumo (nonlinear) system on the axon, coupled with 

• a (linear) Rail model for the dendritical tree, complemented with 

• Kirchhoff-type rule in the soma. 

It is commonly accepted that dendrites conduct electricity in a passive way. The well 
known Rail's model EQl EI] simplify the analysis of this part by considering a simpler, 
concentrated "equivalent cylinder" (of finite length id) that schematizes a dendritical 
tree; he showed that a linear cable equation fits experimental data on dendritical 
trees quite well, provided that it is complemented by a suitable dynamical conditions 
imposed in the interval end corresponding to the soma. Further efforts have been 
put on models for signal propagation along the axon. Shortly after the publication of 
Hodgkin and Huxley's model for the diffusion of electric potential in the squid giant 
axon, a more analytically treatable model was proposed by FitzHugh and Nagumo; the 
model is able to catch the main mathematical properties of excitation and propagation 
using 

o a voltage-like variable having cubic nonlinearity that allows regenerative self- 
excitation via a positive feedback, and 
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o a recovery variable having a linear dynamics that provides a slower negative 
feedback. 

In our model the axon has length £, i.e. the space variable x in the above equations 
ranges in an interval (0,£), where the soma (the cell body) is identified with the point 
0. 

There is a large evidence in the literature that realistic neurobiological models shall 
incorporate stochastic terms to model real inputs. It is classical to model the random 
perturbation with a Wiener process, compare [23] , as it comes from a central limit 
theorem applied to a sequence of independent random variables. 

However, there is a considerable interest in literature to apply different kind of 
noises: we shall mention long-range dependence processes and self-similar processes, 
as their features better model the real inputs: see the contributions in [351 Part II]. 
Further, they can be justified theoretically as they arise in the so called Non Central 
Limit Theorem, see for instance [7J [55] . 

The fractional Brownian motion is of course the most studied process in the class 
of Hermite processes due to its significant importance in modeling. It is not only 
selfsimilar, but also exhibits long-range dependence, i.e., the behaviour of the process 
at time t does depend on the whole history up to time t, stationarity of the increments 
and continuity of trajectories. 

5.1. The abstract formulation 

In the following, as long as we allow for variable coefficients in the diffusion operator, 
we can let the edges of the neuronal network to be described by the interval [0, 1]. The 
general form of the equation we are concerned with can be written as a system in the 
space X = (L 2 (0, l)) 2 x R x L 2 (0, 1) for the unknowns (u, u d , d, v): 



§- t u{t,x) = £(c(x)JLu(t,x)) -p(x)u(t,x)-v(t,x)+d(u(t,x)) + §- t C(t,x) 
■§iU d (t,x) = ^ (c d (x)J^u d (t,x)) -p d (x)u d (t,x) + j^( d {t,x) 
-jd(t) - (c(0)£u(t, 0) - c d (l)JLu d (t, 1)) 



&d(t) 



(5.1) 



§- t v{t,x) =u(t,x)-ev(t,x) + §jC(t,x) 



under the following continuity, boundary and initial conditions 



d(t) = u(t,0) = u d (t,l), 
JL U (t,l) = 0, ^ d (t,0)=0, 
u(0, x) = u (x), v(0,x) = v (x), 



u d (0,x) = u d - (x) 



t > 

i>0 (5.2) 



Throughout the paper we shall assume that the coefficients in (|5.ip satisfy the 
following conditions. 



Assumption 5.1. 

• The function 9 
A > such that 



satisfies some dissipativity conditions: there exists 



for h{u) — — Xu + 9(u) it holds 
[h(u) -h(v)] (u-v) < Vu,wG 



\h(u)\ < c(l + M 2p+1 ), pEN. 



(5.3) 



22 



S. Bonaccorsi, C. A. Tudor 



• c,Cd,p,Pd G C 1 ([0il]) are continuous, positive functions such that, for some 
C>0, 

C<c(x),c d (x) < I C <p(x)-X,p d (x) < ^; 

• 7 > 0, e > are given constants. 

Remark 5.1. The function 9 : R — ► R, in the classical model of FitzHugh, is given 
by 6{u) = u(l - u)(u - £) for some £ G (0, 1); it satisfies with A = |(£ 2 - £ + !)■ 
Other examples of nonlinear conditions are known in the literature, see for instance [5] 
and the references therein. 



Our aim is to write equation (|5.ip . endowed with the conditions in (|5.2p , in an 
abstract form in the Hilbert space X = (L 2 (0, l)) 2 x R x L 2 (0, 1). We also introduce the 
Banach space Y = (C([0, l])) 2 x R x L 2 (0, 1) that is continuously (but not compactly) 
embedded in X. In this section we establish the basic framework that we need in order 
to solve the abstract problem. To this aim we need to prove that the linear part of 
the system defines a linear, unbounded operator A that generates on X an analytic 
semigroup. We shall also study the dissipativity of A and of the nonlinear term F (see 
(El). 



On the domain 
D(A) :-- 



t):=(«,D,(i,iij) T e(fl 2 (0,l)) 2 xIxL 2 (0 I l) s. th. u(0) = u d (l) = d, 

tt / (l)=0, u' d (0)=0, c(0y(0)+c d (iy d (l) = 0J 

(5.4) 



we define the operator A by setting 



Ac := 



f {cu'Y — pu + Au — v \ 

(c d u' d y - p d u d 

- 7 d-( c (0)u'(0)-c d (lK(l)) 
V u — ev J 



(5.5) 



In order to treat the nonlinearity in our system, we introduce the Nemitsky operator 
6 on L 2 (0, 1) such that Q(u)(x) = h(u(x)) for all u G C*([0, 1]) C L 2 (0, 1). Then we 
define F on X by setting 



F(o) = (9( M ),0,0,0) T 
on the domain D(F) = { (u, v, d, u d ) T G X : u G C([0,1])} 

Remark 5.2. In the above setting, the function F satisfies the conditions in Assump- 
tion o 

Finally, setting B(t) = (C(t), ( v (t), 0, ( d {t)) T , we obtain that the initial value 
problem associated with (|5.ip ~ (|5.2[) can be equivalently formulated as an abstract 
stochastic Cauchy problem 

f do(t) = [Ao(*)+F(t>(i))d* + dB(*), t>0, . . 

,(0)^o, (5J) 
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where the initial value is given by Do := (uq, vq, uq(0), u d[ o) T G X. 

In the next section we shall prove that the leading operator A in Eq. (|5.7p satisfies 
the condition in Assumption 11.11 According to Theorem 11.61 this implies that there 
exists a unique solution to problem (|5.7p whenever the noise (B(t))t>o is a fractional 
Brownian motion with Hurst parameter H > i, or a bifractional Brownian motion 
with H > h and K > 1/2H, or an Hermite process with selfsimilarity order H > |, 
or, more generally, a process that satisfies Assumption 11.41 

Theorem 5.2. The proposed model for a neuron cell, endowed with a stochastic input 
that satisfies the conditions in Assumption \l-4\ has a unique solution on the time 
interval [0, T], for arbitrary T > 0, which belongs to 

L%(Q; C([0, T];X)) n L%(fl; L 2 ([0, T]; V)) 

and depends continuously on the initial condition. 
5.2. The well-posedness of the linear system 

As stated above, we can refer to some results in the existing literature in order 
to prove well-posedness and further qualitative properties of our system: the main 
references here are [H [16j US] ■ 

Our first remark is that, neglecting the recovery variable v, the (linear part of the) 
system for the unknown (u, u^, d) is a diffusion equation on a network with dynamical 
boundary conditions: 

Jju(f, x) = ^ (J^c(x)u(t, x)) - p(x)u(t, x) + Xu(t, x) 
^u d (t,x) = ^ (J^c d (x)u d (t,x)) - p d (x)u d (t,x) (5.8) 
&d(t) = - 7 d(t) - (c(O)^tt(t, 0) - c d (l)^u d (t, 1)) 

Such systems are already present in the literature. Let us define X = (L 2 (0, l)) 2 x M 
and introduce the operator 

(cu'Y — pu + Xu 

{c d u' d y -p d u d 

- ll d-(c(0)u'(0)-c d (l)u' d (l)) / 

with coupled domain 

D(A) = {(u, u d , d) T e {H 2 (0, l)) 2 x C : u(0) = u d {l) = d} 
Then, by quoting for instance the papers [TB1[T5|, we can state the following result. 




Proposition 5.3. The operator (A, D(A)) is self-adjoint and dissipative and it has 
compact resolvent; by the spectral theorem, it generates a strongly continuous, analytic 
and compact semigroup (S(t)) t >o on the Hilbert space X . 

The next step is to introduce the operator A on the space X = X x L 2 (0, 1). We 
can think A as a matrix operator in the form 

A -Pi 
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where Pi is the immersion on the first coordinate of X: P\V = (v, 0,0) T , while 
Pi (u,u d ,v) T = u. 

In order to prove the generation property of the operator A, we introduce the Hilbert 
space 



Proposition 5.4. The operator A generates a strongly continuous, analytic semigroup 
(S(t)) t >Q on the Hilbert space X that is uniformly exponentially stable: there exist 
M > I and uo > such that \\S(t)\\ L{x) < Me _ui for all t > 0. 

Proof. We first notice that (A, -D(A)) is the operator associated with the form (a, V): 
compare for instance [H Lemma 4.2]. 

The form (a,V) is non-symmetric, as it can be seen by setting = (1,1,1,0) T 
and u (2) = (1,1,1,1) T and computing a(u (1) ,u (2) ) - a(u (2) ,u (1) ) = 2. However (o,V) 
is densely defined, coercive and continuous, see [H Theorem 4.3]. Then, the properties 
of the semigroup follow from standard theory of Dirichlct forms, compare 19 . 

Notice that the operator A is not self-adjoint, as the corresponding form a is not 
symmetric; also, since V is not compactly embedded in X, it is easily seen that the 
semigroup generated by A is not compact hence it is not Hilbert-Schmidt. For our 
purposes, they are of fundamental importance the following observations. 

Remark 5.3. The form domain V is isometric to the fractional domain power D{{— A) 1 / 2 ). 
This follows since the numerical range of the form a is contained in a parabola, compare 
[H Corollary 6.2], and then by an application of a known result of Mcintosh [HI 
Theorems A and C]. 

Remark 5.4. The form a is real-valued and coercive, hence 



for some to > 0. 

Although we shall not use directly the next result in this paper, we can characterize 
further the specturm of A in the complex plane. This result was first investigated in 
[3]; we provide here our proof for the convenience of the reader. 

Lemma 5.5. The spectrum of A in the complex plane is contained in the union of the 
(discrete, real and negative) spectrum of A and a bounded B. 




c := (u,u d ,d,v) T e (^((Xl)) 2 x R x L 2 (0, 1) s. th. 

u(0) = u d (l) = d 



and the sesquilinear form a : V x V — > K defined by 




(—Au, u) = a(u, u) > oj\\u\\y 



Dissipative stochastic evolution equations driven by general noise 
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Proof. To compute the spectrum we apply [171 Theorem 2.4]. There it is proved 
that for any A ^ a (A) U {-e} it holds A S a (A) if and only if G a(A\(A)), where 
A\(A) is the operator 



(A - A) + 



e + X 



PiP{ 



By standard results on additive bounded perturbations of operators, we notice that 



{A : e a(A x (A))} c D {X 



1 



R(X,A)-—P 1 P^ 



A + e 



< 1} D {A 



1 



1 



\X + lu\ |A + e| 



<1} 



where — u> — s(A) is the spectral bound of A, that is a negative real number by 
Proposition 15.31 Therefore, setting B = {X : \X + ui\ \X + e\ < 1} we have that B 
is a bounded subset of the complex plane. 
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